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1.  SUMMARY 

The  development  of  a  high-speed,  phase- locked,  real¬ 
time,  point  diffraction  interferometry  system  for  quan¬ 
titative  imaging  unsteady  separated  flows  is  described. 
The  system  enables  recording  of  up  to  224  interfer- 
ograms  of  the  dynamic  stall  flow  over  an  oscillating 
airfoil  using  a  drum  camera  at  rates  of  up  to  40  KHz 
controlled  by  custom  designed  electronic  interlocking 
circuitry.  Several  thousand  interferograms  of  the  flow 
have  been  obtained  using  this  system.  A  comprehen¬ 
sive  image  analysis  package  has  been  developed  for 
automatic  processing  of  this  large  number  of  images. 
The  software  has  been  specifically  tuned  to  address  the 
special  characteristics  of  airfoil  flow  interferograms. 
Examples  of  images  obtained  using  the  standard  and 
the  high-speed  interferometry  techniques  are  presented 
along  with  a  demonstration  of  the  image  processing 
routine’s  ability  to  resolve  the  fine  details  present  in 
these  images. 

LIST  OF  SYMBOLS 

Cp  pressure  coefficient 

c  airfoil  chord 

f  frequency  of  oscillation,  Hz 

k  reduced  frequency  = 

I  image  intensity 

M  freest  ream  Mach  number 

Uoo  freest  ream  velocity 

x,y  chordwise  and  vertical  distance 

a  angle  of  attack 

ao  mean  angle  of  attack 

Oi  amplitude  of  oscillation 

€  fringe  number 

p  density 

Pq  density  ‘at  atmospheric  conditions 

pr  density  at  reference  conditions 


2.  INTRODUCTION 

Dynamic  stall  has  limited  the  flight  envelope  of  heli¬ 
copters  for  many  years.  The  problem  has  been  stud¬ 
ied  both  in  the  laboratory  and  in  flight,  employing 
complex  surface  measurement  techniques  such  as  pres¬ 
sure  transducers  or  skin  friction  gauges.  To  properly 
understand  the  complex  physics  associated  with  dy¬ 
namic  stall  both  surface  and  off-surface  measurements 
are  needed  simultaneously.  Quantitative  visualization 
of  the  flow  field  during  compressible  conditions  can  be 
used  in  such  cases,  but  in  the  past  this  heis  been  pos¬ 
sible  only  through  carefully  aligned  and  meticulously 
reconstructed  holographic  interferometry. 

As  part  of  a  long-range  effort  focused  on  exploring 
the  physics  of  compressible  dynamic  stall,  a  research 
wind  tunnel  was  developed  at  NASA  Ames  Research 
Center(ARC)  which  permits  visual  access  to  the  full 
flow  field  surrounding  an  oscillating  airfoil  during  com¬ 
pressible  dynamic  stall  (Ref.  1).  Initially,  a  strobo¬ 
scopic  schlieren  technique  was  used  for  visualization 
of  the  stall  process  (Ref.  2),  but  the  primary  research 
tool  has  been  the  real-time,  self-aligning  technique  of 
Point  Diffraction  Interferometry(PDI),  which  has  been 
carefully  optimized  for  use  in  this  project  (Ref.  2,3). 
PDI  is  a  robust,  self- aligning  interferometry  technique. 
It  uses  conventional  schlieren  optical  components  and 
is  immune  to  the  problems  that  generally  limit  the  ap¬ 
plication  of  other  interferometers.  Using  PDI,  both 
global  and  surface  flow  features  can  be  simultaneously 
documented.  It  has  been  shown  in  Ref.  2  that  dy¬ 
namic  stall  onset  occurs  over  a  very  narrow  angle  of 
attack  range  (0(0.5)  deg),  at  instantaneous  angles  of 
attack  when  the  airfoil  pitch  rate  is  mciximum,  dur¬ 
ing  which  a  large  leading  edge  vortex  forms  and  grows 
rapidly  from  a  fully  attached  flow.  Further,  hysteresis 
results  in  cycle-to-cycle  variation  of  stall  onset  which 
smears  the  data  if  acquired  by  conventional  ensemble 
averaging.  Thus,  it  is  required  that  the  flow  devel¬ 
opments  be  captured  in  just  one  pitch-up  cycle  which 
places  extraordinary  demands  on  the  data  acquisition 


^  Mailing  Address:  M.S.  260-1,  NASA  Ames  Research  Center,  Moffett  Field,  CA  94035,  U.S.A. 


21-1 


* 


methods.  Since  the  pitch  rate  could  be  as  large  as  4000 
deg/sec  and  the  data  is  needed  at  a  very  high  resolu¬ 
tion  to  document  the  details  of  dynamic  stall  evolu¬ 
tion,  extremely  rapid  real-time  imaging  rates  {O  f20 
KHz))  become  necessary.  It  is  clear  that  holographic 
techniques  are  not  applicable  here. 

The  recording  of  the  rapid  flow  changes  also  requires 
freezing  the  flow  events  during  data  acquisition.  This 
implies  sampling  times  of  less  than  a  microsecond  for 
the  present  study.  The  energy  output  of  even  the  most 
powerful  laboratory  laser  is  generally  less  than  200 
for  repetitively  pulsed  conditions.  Imaging  at  such  a 
low  light  level  restricts  the  choice  of  recording  sys¬ 
tems  to  photographic  film  since  even  high  speed  videos 
presently  do  not  have  a  satisfactory  response  at  these 
low  light  levels.  Trials  showed  that  such  systems  re¬ 
sponded  randomly  to  individual  photons  of  light  rather 
than  to  light  across  the  whole  field  of  view. 

One  of  the  most  valuable  aspects  of  PDI  is  the  fact 
that  interferograms  can  be  produced  in  real-time  on  a 
continuous  basis  through  the  use  of  a  rapidly-pulsed 
laser.  Pulsing  at  rates  of  up  to  40  KHz  and  record¬ 
ing  with  high-speed  film  camera  produces  interfero¬ 
grams  at  a  sufficient  resolution  (both  in  space  and 
time)  to  analyse  the  rapidly  developing  unsteady  dy¬ 
namic  stall  field.  In  the  present  experiments,  this  has 
been  achieved  through  the  use  of  a  high-speed  drum 
camera  combined  with  custom  designed  electronic  cir¬ 
cuitry.  Interferogram  sequences  of  dynamic  stall  devel¬ 
opment  during  a  single  oscillation  cycle  of  the  airfoil 
have  been  recorded  at  rates  of  up  to  20  KHz.  A  de¬ 
tailed  analysis  of  the  effects  of  pitch  rate,  Mach  num¬ 
ber,  Reynolds  number,  amplitude  of  oscillation,  and 
other  parameters  on  the  dynamic  stall  process  has  thus 
become  possible. 

The  present  research  has  also  focused  on  quantitative 
determination  of  the  fluid  physics  of  the  compress¬ 
ible  dynamic  stall  flow  field,  including  the  gradients  of 
pressure  in  space  and  time  from  these  interferograms. 
The  development  of  image  analysis  techniques  specific 
to  PDI  interferograms  obtained  here  is  a  major  effort 
in  this  quantitative  evaluation.  Instantaneous  pressure 
distributions  can  now  be  obtained  semi-automatically 
making  practical  the  analysis  of  the  thousands  of  in¬ 
terferograms  that  are  produced  in  this  research. 

3.  DESCRIPTION  OF  THE  FACILITY  AND 
MEASUREMENT  TECHNIQUE 

3.1.  The  Compressible  Dynamic  Stall  Facility 

The  experimental  studies  were  carried  out  in  the  Com¬ 
pressible  Dynamic  Stall  Facility (CDSF)  of  the  ARC 
Fluid  Mechanics  Laboratory.  The  CDSF  is  an  indraft 
wind  tunnel  with  a  10  in  x  14  in  test  section  and  is 
equipped  with  a  drive  for  producing  a  sinusoidal  vari¬ 
ation  of  the  airfoil  angle  of  attack.  The  tunnel  Mach 
number  is  controlled  by  a  choked,  variable- area,  down¬ 
stream  throat  in  the  range  of  0.1  <  M  <  0.5.  The  flow 
is  produced  by  a  6MW,  240,000  CFM  continuously 
running  evacuation  compressor.  The  airfoil  mean  an¬ 
gle  of  attack  can  be  set  to  0  <  a  <  15°,  the  amplitude 
of  oscillation  to  2°  <  ai  <  10°,  and  the  oscillation 
frequency  to  0  <  /  ^  lOOHz,  In  the  CDSF  a  3-inch 


chord  NACA  0012  airfoil  is  supported  by  pins  between 
two  6  in  diameter  optical  glass  windows.  This  unique 
mounting  method  provides  direct  optical  access  to  the 
airfoil  surface  permitting  flow  exploration  using  nonin- 
trusive  diagnostic  techniques,  Fig.  1.  Tests  on  a  6-inch 
chord  airfoil  mounted  between  metal  ports,  with  glass 
inserts  for  optical  access  to  the  first  40%  chord  of  the 
airfoil  are  used  to  establish  Reynolds  number  effects 
produced  at  a  given  freestream  Mach  number.  Block¬ 
age  effects  due  to  the  larger  airfoil  are  small  around 
dynamic  stall  onset  angles  of  attack.  (The  correction 
to  pressures  is  less  than  5%  at  a  =  10  deg.)  Encoders 
on  the  drive  system  record  the  mean  and  instantaneous 
angles  of  attack.  The  encoder  signals  are  processed  in 
a  custom  built  interface  known  as  the  Oscillating  Air¬ 
foil  Position  Interface(OAPI)  which  is  also  used  for 
conditional  sampling  and  phase  locking. 

3.2.  Description  of  the  PDI  Technique 

3. 2.  A.  The  Interferometry  System 

As  stated  earlier,  the  experimental  studies  used  the 
real-time  technique  of  point  diffraction  interferometry. 
PDI  provides  detailed,  instantaneous  and  quantitative 
flow-field-density  information,  from  which  both  sur¬ 
face  and  global  pressure  distributions  can  be  derived. 
The  PDI  optical  arrangement  is  similar  to  that  of  a 
schlieren  system,  but,  as  shown  in  Fig.  2,  a  laser  light 
source  is  used  with  a  beam  expander  to  fill  the  entire 
field  of  view  (determined  by  the  tunnel  windows)  and 
a  point  diffractor  (pin-hole)  replaces  the  knife  edge. 
The  optics  are  aligned  in  the  standard  Z-type  config¬ 
uration  to  minimize  astigmatism.  Imaging  optics  are 
set  up  further  downstream  along  the  beam  path  for 
recording  the  flow. 

The  technique  uses  one  single  pass  of  the  laser  beam 
through  the  test  section  and  depends  upon  the  ability 
of  a  pin-hole  created  in-siiu  in  a  semi-transparent  plate 
to  produce  the  reference  beam.  Creation  of  the  pin¬ 
hole  requires  focusing  the  laser  beam  on  the  plate  and 
exposing  the  plate  emulsion  to  a  very  high  level  of  laser 
energy  with  no-flow  in  the  tunnel.  As  the  photographic 
emulsion  burns  away,  a  clear  spot  (pin-hole)  appears 
on  the  plate.  With  the  flow-on,  light  phase  shifted  by 
the  flow  density  changes  (signal  beam)  passes  around 
this  spot  to  produce  interference  fringes  on  a  continu¬ 
ous  basis  in  real-time.  The  portion  of  the  beam  passing 
through  the  pin-hole  becomes  the  reference  beam  due 
to  the  spatial  filtering  characteristics  of  the  pin-hole. 
The  technique  has  now  evolved  considerably  and  sev¬ 
eral  hundred  interferograms  can  be  obtained  in  a  day’s 
work. 

3.2. B.  Details  of  the  High-Speed  Camera,  Laser  Con¬ 
trol  and  Recording  Technique 

A  Qunatronix  Series  100  CW/pumped  Nd:YAG  laser, 
capable  of  operating  from  DC  to  50  KHz  wais  used  in 
the  experiments.  It  could  be  externally  triggered  with¬ 
out  any  detectable  delay  at  all  rates.  The  pulse  dura¬ 
tion  and  the  energy  output  varied  nonlinearly  from  85 
ns  and  140  //J  at  500  Hz,  420  ns  and  25  /iJ  at  40  KHz 
and  100  ns  and  11  at  50  KHz.  At  the  rates  used  for 
the  high-speed  interferometry  experiments  being  re- 
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ported,  the  corresponding  numbers  were:  140  ns  and 
65  /iJ  at  10  KHz  and  240  ns  and  17  at  20  KHz,  at 
nearly  full  current  settings.  The  energy  density  in  the 
laser  light  pulse  at  the  10  KHz  rate  was  adequate  to 
give  proper  exposure  on  ASA  100  T-MAX  film;  ASA 
400  film  was  necessary  at  20  KHz. 

A  variable  speed  Cordin  drum  camera  (DYNAFAX 
Model  350)  was  used  for  image  recording.  A  rotating 
8-faceted  mirror  in  the  camera  reflected  the  incoming 
light  beam  onto  film  rotating  in  the  same  direction  in 
the  camera  drum.  Effective  shutter  times  of  1.35  fisec 
could  be  achieved  at  the  40  KHz  framing  rate.  At  20 
KHz,  this  time  was  2.7  fisec.  The  camera  recorded 
two  rows  of  16  mm  images  on  a  35  mm  film  strip,  with 
successive  exposures  recorded  alternately  in  each  row, 
but  displaced  by  16  frames.  A  maximum  of  224  frames 
could  be  recorded  at  any  framing  speed. 

The  laser  was  triggered  by  TTL  pulses  emitted  by  cus¬ 
tom  designed  and  built  (in-house)  circuitry  installed 
on  the  camera.  An  infrared  (IR)  emitter/detector  was 
installed  in  the  camera  (out  of  the  light  path  as  shown 
in  Fig.  3)  to  detect  reflections  from  the  mirror  facets 
as  they  passed.  The  selected  film  was  not  sensitive  to 
the  940nm  IR  wavelength.  Each  mirror  facet  detection 
produced  two  pulses  using  fast  rise-time  (1  nsec)  photo 
diodes  to  ensure  adequate  signal  level.  The  pulses  were 
delayed  from  the  time  of  the  IR  detection  pulse.  The 
time  delays  were  calibrated  so  that  the  mirror  facet 
was  aligned  with  the  image  frame  at  the  time  of  each 
laser  pulse.  Two  photo  diodes  were  placed  in  the  cam¬ 
era,  one  at  the  frame  position  in  each  film  track  for  cal¬ 
ibrating  the  time  delays,  then  removed  for  operation. 
The  tuning  procedure  involved  adjusting  two  delay 
times  with  the  camera  running:  Ti ,  the  delay  between 
detecting  a  mirror  facet  and  emitting  the  first  TTL 
pulse  (the  trigger  pulse  to  the  laser  and  for  data  col¬ 
lection)  and  T2,  the  time  between  the  two  TTL  pulses 
(the  time  between  frames).  The  delay  times  Ti  and  T2 
were  adjusted  to  maximize  the  laser  light  detected  by 
the  frame  photo  detectors.  Once  tuning  was  properly 
completed,  the  photo  detectors  were  moved  from  the 
field  of  view  to  permit  laser  light  to  reach  the  film 
plane.  The  short  effective  shutter  times  (of  1.35  /isec 
at  the  maximum  camera  speed)  and  the  high  framing 
speeds  required  a  careful  design  of  the  electronic  sys¬ 
tem  that  included  schemes  for  proper  attenuation  of 
noise. 

The  lenses  in  the  PDI  imaging  optics  could  be  adjusted 
to  produce  a  focused  image  of  the  flow  area  at  the  film 
plane  of  the  camera.  Aligning  the  camera  along  the 
optical  axis  of  the  interferometry  system  required  very 
accurate  adjustment. 

3.2. C.  Operation 

The  interframe  pulse  delay  was  tuned  to  the  desired 
rate  and  the  actual  rate  of  the  camera  was  measured 
using  a  frequency  counter.  Several  interferogram  se¬ 
quences  were  then  obtained  at  11.56  KHz  and  19.62 
KHz  for  varying  initial  phase  angles.  In  order  to  main¬ 
tain  a  consistent  pulse  energy  level,  the  laser  was  trig¬ 
gered  by  an  external  pulse  train  at  a  40  KHz  rate  be¬ 
fore  the  images  were  acquired.  This  was  necessary 
to  protect  the  laser  crystal  from  the  giant  pulse  that 


is  normally  generated  when  the  laser  is  pulsed  after  a 
short  lapse  time.  In  order  to  prevent  these  pulses  from 
exposing  the  film  and  for  safety  reasons,  a  solenoid  ac¬ 
tuated  laser  shutter  was  set  up  in  front  of  the  laser.  A 
hand  switch  Wcis  used  to  initiate  the  controlled  laser 
pulsing  sequence,  which  is  schematically  described  in 
Fig.  4.  The  corresponding  timing  sequence  is  shown 
in  Fig.  5.  After  the  switch  was  pressed,  the  circuitry 
wcis  activated  by  an  event  pulse  from  the  OAPI  corre¬ 
sponding  to  the  manually  preselected  angle  of  attack 
which  in  turn  triggered  the  laser  safety  shutter.  The 
laser  pulsing  circuit  was  then  inhibited  (for  1.1  msec) 
until  the  laser  shutter  fully  opened.  The  laser  was 
enabled  at  the  expiration  of  the  delay  and  was  actu¬ 
ally  triggered  from  the  next  camera  pulse,  at  which 
time  the  encoder  was  latched  and  recorded  in  a  512 
word  first-in-first-out(FIFO)  buffer.  During  this  short 
elapsed  time,  the  laser  built  up  sufficient  charge  to 
cause  the  first  pulse  to  be  a  ‘‘small”  giant-pulse,  which 
over-exposed  the  first  frame.  This  frame  served  to 
identify  the  first  image  on  the  film  strip;  thus,  it  was 
possible  to  accurately  match  the  interferogram  images 
with  the  phase  angle  of  motion  and  to  correlate  the 
values  in  the  FIFO  buffer.  A  frame  counter,  started  at 
the  first  laser-pulse  event  inhibited  the  laser  after  200 
laser  pulses.  Each  film  strip  can  record  a  maximum 
of  224  images.  Following  the  completion  of  the  imag¬ 
ing,  the  shutter  was  closed  and  the  laser  returned  to 
the  constant  40  KHz  external  triggering.  The  camera 
alignment  was  verified  by  taking  test  sequences  on  a 
Polapan  ASA  125  film;  the  data  was  obtained  on  the 
higher  resolution  T-MAX  400  film. 

3.2.D.  Determination  of  Pressure  Coefficients  from 
Fringe  Number 

The  quantitative  nature  of  the  interferograms  allowed 
derivation  of  the  pressure  distribution  over  the  airfoil 
when  the  flow  was  attached.  The  interferograms  were 
processed  in  a  manual  mode  using  a  specially  devel¬ 
oped  software  package.  The  program  read  a  digitized 
(256  gray  level)  interferogram  on  an  IRIS  Work  Station 
and  over-layed  an  airfoil  using  the  registration  markers 
on  the  images.  The  intersection  of  the  fringes  with  the 
airfoil  upper  and  lower  surfaces  (or  the  local  boundary 
layer  edge,  when  detectable)  were  interactively  picked 
by  the  user.  Since  each  fringe  is  a  line  of  constant 
density,  the  corresponding  pressure  at  the  boundary 
layer  edge  could  be  calculated  using  isentropic  flow  re¬ 
lations.  This  pressure  was  then  used  as  the  surface 
pressure,  under  the  boundary  layer  assumptions.  The 
density  along  any  fringe  could  be  calculated  from  the 
Gladstone  -  Dale  equation,  which  for  the  present  wind 
tunnel  and  laser  simplifies  to 


p-pr  =  0.009421e 


where  e,  the  fringe  number,  is  0,i:l,±2,...  for  the 
bright  fringes  and  for  the  dark  fringes. 

Fringes  from  the  free  stream  to  the  stagnation  point 
have  positive  values.  The  corresponding  Cp  values 
were  then  computed  from  the  relation 


21-3 


- 1' 

=  — H- 

In  the  end,  the  program  provided  an  output  data  file 
containing  the  various  physical  variables,  in  a  format 
suitable  for  plotting.  Typical  processing  time  was 
about  3-5  minutes  per  image.  In  cases  where  the 
fringe  density  was  high  or  the  fringes  were  fuzzy,  the 
user  could  go  into  the  ‘off-body’  mode  and  pick  fringes 
along  a  line  parallel  to  and  away  from  the  airfoil  sur¬ 
face  where  the  fringes  are  farther  apart.  For  this  pur¬ 
pose,  an  option  to  superpose  two  larger  airfoils  over  the 
image  on  the  screen  was  provided.  The  fringe  intersec¬ 
tions  on  the  larger  airfoils  were  then  suitably  projected 
on  to  the  airfoil  surface.  At  angles  of  attack  near  the 
dynamic  stall  angle,  the  fringes  near  the  leading  edge 
region  were  very  dense  reflecting  the  large  local  density 
gradients.  Further,  in  this  region,  optical  noise  intro¬ 
duced  by  the  shadowgraph  effect  generally  lowered  the 
contrast,  making  it  a  location  where  the  off-body  mode 
needed  to  be  invoked. 

In  the  present  study  the  entropy  change  in  the  vor¬ 
tical  flow  was  ignored  (for  lack  of  a  better  method). 
Interferograms  with  shocks  have  not  been  processed 
for  pressures  because  of  this  limitation. 

3.3.  Details  of  Image  Processing 

Aside  from  these  manual  procedures,  a  program  was 
also  developed  for  automatic  analysis  of  the  interfer¬ 
ograms.  A  number  of  digital  filters  and  image  en¬ 
hancement  procedures  are  available.  The  operator  is 
required  to  register  the  image  (thus  masking  the  air¬ 
foil  and  establishing  a  coordinate  system),  choose  the 
filtering  and/or  enhancement  methods,  and  mark  re¬ 
gions  of  the  image  for  special  treatment.  Fringe  cen¬ 
terlines  are  automatically  traced  with  processing  times 
approximately  one  quarter  that  of  fully-manual  pro¬ 
cessing,  and  operator  bias  (a  potential  source  of  error 
in  the  tracing  of  fringes)  is  eliminated. 

Full  automation  of  image  processing  requires  software 
capable  of  addressing  the  special  characteristics  of 
these  images,  such  as  variations  in  lighting  and  con¬ 
trast  across  an  image  that  are  inherent  to  the  tech¬ 
nique,  an  occasional  broken  (or  split)  fringe,  the  high 
fringe  density  near  the  suction  peak  (as  many  as  50 
fringes/mm  along  the  airfoil  surface),  presence  of 
shocks,  etc.  Full  details  of  the  image  processing 
method,  e.g.  image  filtering,  fringe  centerline  detec¬ 
tion,  fringe  identification,  etc.,  can  be  found  in  Ref.  5. 
The  flow  chart  shown  in  Fig.  6  illustrates  the  typical 
processing  procedure.  The  need  for  the  spin  filter,  an 
intensity  gradient  aligned  spatial  filter  (Ref.  7),  de¬ 
pends  on  the  quality  of  the  original  image.  Filtering 
and/or  contrast  enhancement  are  perfornied  primarily 
to  improve  edge  detection,  the  first  stage  in  fringe  cen¬ 
terline  detection,  as  described  below.  For  the  results 
presented  here,  the  region  around  the  leading  edge  of 
the  airfoil  was  magnified  and  treated  separately  from 
the  rest  of  the  image. 

Fringe  centerlines  are  located  by  identifying  the  local 
extrema  in  the  image  intensity.  These  extrema  are 


identified  within  the  background  noise  through  a  two- 
stage  procedure  which  first  locates  the  fringe  edges, 
and  then  fits  a  smooth  polynomial  to  the  data  be¬ 
tween  successive  edges.  The  location  of  the  extreme 
value  of  the  approximated  data  is  taken  as  the  loca¬ 
tion  of  the  fringe  centerline.  This  procedure  is  a  one 
dimensional  operation  performed  in  two  passes  over 
the  image;  rows  of  pixels  are  examined  in  the  first  pass 
and  columns  are  examined  in  the  second  pass. 

A  zero-crossing  algorithm  is  used  to  identify  the  fringe 
edges.  The  zero-  crossing  threshold  level  is  either  the 
mean  intensity  along  the  row  or  column  of  pixels  un¬ 
der  consideration,  or  the  global  mean  intensity  of  the 
image.  The  global  mean  is  used  if  the  current  row  or 
column  follows  a  fringe  or  crosses  only  one  edge.  A 
noise  band  around  the  threshold  is  specified  to  elimi¬ 
nate  hysteresis  crossings  due  to  noise  in  the  image. 

The  edges  are  used  as  a  guide  to  locating  the  fringe 
centerlines.  Assuming  that  each  successive  pair  of 
edge  points  outline  a  single  fringe,  the  centerline  of  the 
fringe  will  pass  through  the  point  of  extreme  intensity 
between  these  two  edge  points.  In  order  to  unambigu¬ 
ously  locate  this  extreme,  a  smooth  polynomial  curve 
is  fit  to  the  data  between  pairs  of  edge  points  using 
Chebyshev  approximation.  The  image  intensity  is  ap¬ 
proximated  on  the  interval  [-1,  1]  as 

m  —  1  I 

k  =  0 

where  the  Chebyshev  polynomials  and  the  coefficients 
are  given  by 

Tk{x)  =  cos{k  cos^^x) 

and 

jzzl 

respectively,  and  the  data  {I{x))  and  the  Chebyshev 
polynomials  (Tjb(a:))  are  evaluated  at  the  N  zeros  of 
Tn{x),  namely 


i  =  . ''■ 

The  coefficients  are  calculated  for  Chebyshev  polyno¬ 
mials  of  degree  N,  while  the  approximating  polynomial 
is  truncated  to  a  polynomial  of  degree  m  ^  N  (typi¬ 
cally  AT  =  20,  m  =  3).  This  yields  an  accurate  approx¬ 
imation  of  degree  m  to  the  data  [8j.  For  sufficiently 
large  N,  asm N  the  approximating  polynomial  re¬ 
produces  the  original  data,  including  the  noise.  Due 
to  the  multiple  inflection  points  in  the  data,  approx¬ 
imating  data  spanning  more  than  one  fringe  requires 
a  higher  order  polynomial  (m),  which  tends  to  repro¬ 
duce  the  noise  content.  Therefore,  the  approximation 
procedure  is  confined  to  a  single  fringe. 

This  procedure  was  tested  on  simulated  fringe  data  (a 
modulated  cosine  function  with  added  Gaussian  noise) 
and  was  found  to  be  accurate  to  within  2  pixels  for 


21-4 


fringes  from  5  to  40  pixels  in  width  and  to  within  5 
pixels  for  fringes  up  to  240  pixels  wide. 

The  points  located  this  way  comprise  two  sets  of  pix- 
els,  one  associated  with  the  dark  fringe  centerlines  and 
the  other  with  the  light  fringe  centerlines.  These  data 
sets  form  two  search  lists  which  are  sorted  into  smaller 
groups,  each  tracing  a  single  fringe.  Tracing  begins  ar¬ 
bitrarily  with  the  first  point  in  the  search  list.  A  near¬ 
est  neighbor  search  is  performed  with  an  inner  and 
outer  search  radius.  Any  point  within  the  inner  search 
radius  of  the  current  point  is  automatically  designated 
as  belonging  to  the  current  fringe.  Any  point  found  be¬ 
tween  the  inner  and  outer  search  radii  are  assigned  to 
the  current  fringe  if  certain  criteria  are  satisfied.  These 
criteria  are  designed  to  prevent  the  tracing  procedure 
from  bridging  adjacent  fringes  by  examining  the  im¬ 
age  intensity  across  the  gap  and  by  checking  for  the 
presence  in  the  gap  of  data  points  belonging  to  the 
opposite  fringe  type.  Each  data  point  assigned  to  a 
fringe  is  removed  from  the  search  list.  The  search  is 
then  repeated  relative  to  each  new  point  added  to  the 
fringe  in  the  previous  search.  The  search  ends  when 
the  end  of  the  fringe  is  reached  or  when  a  gap  larger 
than  the  outer  search  radius  is  encountered.  A  new 
fringe  is  designated  and  the  tracing  procedure  begins 
anew  with  the  first  point  in  the  search  list  becoming 
the  first  point  in  the  new  fringe.  The  tracing  is  pro¬ 
gressively  fcLster  2LS  data  points  are  removed  from  the 
search  list  and  placed  in  fringe  lists. 

The  inner  search  radius  is  two  pixels,  hence,  any  one- 
pixel  gaps  in  the  data  are  automatically  bridged.  The 
outer  radius  varies  with  fringe  spacing  and  is  six  pix¬ 
els  in  the  outer  flow  regions  where  the  fringe  spacing 
is  wide,  and  reduces  to  three  pixels  near  the  leading 
edge  where  fringe  spacing  is  often  less  than  five  pixels. 
Due  to  the  small  fringe  spacing  in  this  region  any  noise 
in  the  data  can  cause  the  fringe  tracing  algorithm  to 
erroneously  join  fringes  together.  Therefore,  the  lead¬ 
ing  edge  is  generally  treated  separately  in  the  form  of 
a  magnified  image  segment. 

3.4.  Experimental  Uncertainties 

The  uncertainty  in  Cp  depends  on  the  fringe  number 
under  consideration  and  is  1  fringe  for  the  flow  in  gen¬ 
eral  with  about  3  fringes  possibly  undetectable  near 
the  suction  peak  at  M  =  0.3.  Since  correction  for  solid 
and  wake  blockage  was  less  than  5%  for  Cp  =  —6.0 
at  M  —  0.3,  only  uncorrected  PDI  derived  pressures 
are  reported.  The  losses  in  the  tunnel  screens  caus¬ 
ing  a  decre2ise  in  the  stagnation  pressure  have  been 
included  in  the  computation  of  the  reference  density 
in  this  otherwise  atmospheric  flow  wind  tunnel. 

By  far  the  largest  uncertainty  appears  from  the  very 
strong  density  gradients  near  the  airfoil  leading  edge 
which  introduces  strong  refractive  index  gradients. 
The  effect  of  this  is  to  deflect  the  light  rays  towards 
higher  density  regions.  Hence,  a  dark  region  develops 
on  the  upper  surface  near  the  leading  edge  (Ref.  6) 
in  some  images.  No  corrections  were  applied  to  the 
Cp  values  on  this  account  because  of  the  strong  de¬ 
pendence  of  the  local  density  on  instantaneous  angle 
of  attack,  freestream  Mach  number,  and  reduced  fre¬ 
quency.  Also,  it  was  found  that  by  slightly  realigning 


the  PDI  optics,  much  of  the  distortion  could  be  re¬ 
duced.  With  this  solution,  it  was  estimated  that  the 
true  suction  peak  location  could  be  displaced  by  less 
than  0.5%  of  airfoil  chord  towards  the  trailing  edge. 
Support  for  this  can  be  found  by  directly  comparing 
the  pressures  meaisured  using  surface  taps  and  those 
obtained  with  PDI  as  shown  in  Fig.  7  (Ref.  9)  for 
the  6-inch  airfoil  in  steady  flow  at  M  =  0.3.  At  a  =  8 
deg,  Fig.  7a,  the  two  agree  very  well.  At  a  =  13  deg, 
Fig.  7b,  there  is  agreement  to  within  the  one  fringe 
standard  uncertainty  of  the  PDI  method.  The  differ¬ 
ences  seen  at  the  laminar  separation  bubble  location 
are  due  to  the  insufficient  number  pressure  taps  over 
the  airfoil  that  make  it  impossible  to  record  the  pres¬ 
sures  at  the  resolution  of  PDI. 

With  these  considerations  in  mind,  the  following  are 
the  estimated  uncertainties  in  the  data: 

Mach  number:  ±0.005 

angle  of  attack:  0.05  degrees 

reduced  frequency:  0.005 


-0.45  at  M  =  0.4 

4.  RESULTS  AND  DISCUSSION 

Figure  8  presents  two  interferograms  obtained  in  the 
studies  at  M  =  0.3,  k  =  0.05  and  o  =  10  deg.  Fig. 
8a  shows  a  standard  PDI  interferogram,  whereas  that 
in  Fig.  8b  was  obtained  using  the  high-speed  PDI 
imaging  method  at  a  filming  rate  of  11.56  Khz.  The 
sharpness  of  the  image  in  Fig.  8a  indicates  the  quality 
of  the  pictures  that  can  be  obtained  with  PDI.  The 
high-speed  image  in  Fig.  8b  is  still  clean  even  after  a 
linear  magnification  of  about  35.  The  figure  shows  a 
stagnation  point  on  the  lower  surface  and  as  the  flow 
accelerates,  a  laminar  separation  bubble  forms  down¬ 
stream  of  the  suction  peak.  The  separation  bubble 
can  be  inferred  where  the  fringes  running  parallel  to 
the  upper  surface  abruptly  turn  normal  to  the  surface. 
The  parallel  fringes  represent  the  plateau  in  the  pres¬ 
sure  distributions  due  to  the  presence  of  the  bubble.  In 
both  figures,  a  slight  dark  bulge  is  noticeable  near  the 
airfoil  leading  edge,  which  is  attributed  to  the  steep 
local  density  gradient.  The  global  pressure  maps  are 
shown  in  Fig.  9  for  the  airfoil  at  an  instantaneous  angle 
of  attack  of  13.99  deg.  The  comparison  also  clearly  il¬ 
lustrates  the  differences  that  can  be  present  in  the  flow 
from  one  oscillation  cycle  to  the  next  at  angles  near  the 
dynamic  stall  onset  angle.  These  differences  are  espe¬ 
cially  critical  when  considering  dynamic  stall  control. 
Our  work  has  shown  that  dynamic  stall  for  these  test 
conditions  arises  from  the  bursting  of  a  laminar  sep¬ 
aration  bubble,  which  occurs  when  the  instantaneous 
pressure  gradient  exceeds  a  critical  value  and  so  the 
envelope  of  activation  of  a  control  devise  for  manip¬ 
ulating  the  adverse  pressure  gradient  will  change  for 
each  cycle. 

Figure  10  shows  the  effects  of  pin-hole  size  on  the  PDI 
images  at  M  =  0.4.  Evaluation  of  the  pin-hole  requires 
recording  PDI  images  on  polaroid  film.  The  photo¬ 
graphic  plate  needs  to  be  pre-developed  to  a  level  of 
darkness  that  results  in  an  optical  density  of  around 
1.0.  The  PDI  images  were  found  to  be  generally  insen- 
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sitive  to  the  shape  and  size  of  the  P DI  spot  in  the  range 
tested  as  long  as  the  spot  is  smaller  than  the  Airy  disk 
diameter  of  the  laser  beam.  Good  interferograms  were 
obtained  with  pin-holes  created  by  images  of  both  the 
nearly  circular  windows  and  the  L-shaped  optical  glass 
inserts.  The  most  noticeable  effect  of  a  large  pin-hole 
is  the  many  split  fringes  that  form  as  seen  in  Fig.  10. 
These  can  lead  to  an  incorrect  counting  of  the  fringes, 
causing  erroneous  pressure  distributions  to  be  derived 
from  the  automated  processes.  It  should  be  noted  that 
satisfactory  interferograms  have  been  obtained  in  these 
studies  over  a  range  of  pin-hole  sizes.  Also,  somewhat 
large  pin-holes  are  needed  as  the  flow  Mach  number 
increases  due  to  the  large  density  gradient  and  hence 
a  larger  deflection  of  the  laser  light  near  the  leading 
edge.  Images  can  be  obtained  reliably  without  much 
effort  if  a  PDI  plate  with  the  right  optical  density  is 
used  for  each  Mach  number.  In  general,  it  has  been 
found  that  a  dark  plate  (optical  density  of  the  0(1)) 
is  needed  to  properly  capture  the  finer  flow  details  for 
the  Mach  numbers  in  the  range  of  0.3-0.45.  Lighter 
plates  worked  well  at  lower  Mach  numbers. 

Figures  11  and  12  show  the  results  that  have  been  ob¬ 
tained  with  the  automated  image  processing  program 
for  M  =  0.3  and  M  =  0.45  respectively.  The  image 
in  Fig.  11a  shows  the  airfoil  at  a  =  8.84  deg  for  M 
=  0.3  and  k  =  0.05.  The  corresponding  map  of  the 
dark  fringes  is  plotted  in  Fig.  11b.  Fig.  12a  and  12b 
show,  respectively,  an  image  and  its  fringe  map  ob¬ 
tained  at  M  =  0.45,  k  =  0.05  and  a  =  10^  on  the 
upstroke.  Under  these  conditions  the  local  velocity  ex¬ 
ceeded  the  sonic  speed  and  an  oblique  shock  formed. 
Moreover,  the  interaction  between  the  shock  and  the 
laminar  boundary  ayer  (Re  =  9  x  10^)  resulted  in  the 
formation  of  multiple  shocks.  An  enlarged  view  of  the 
leading  edge  is  shown  to  emphasize  the  shocks.  The 
image  was  magnified  by  three  and  was  contrast  en¬ 
hanced.  No  filtering  was  required.  Despite  the  tight 
fringe  spacing,  and  distortions  due  to  the  shocks  con¬ 
siderable  success  has  been  achieved  in  processing  the 
image  without  manual  intervention.  Typical  process¬ 
ing  times  were  about  5  min.  with  complicated  images 
requiring  longer  times. 

5.  CONCLUSIONS 

A  novel  high-speed,  real-time,  phase- locked,  interfer¬ 
ometry  system  has  been  developed 
and  applied  to  study  oscillating  airfoil  compressible 
dynamic  stall  flow  fields.  The  technique  can  generate 
very  sharp  images  of  the  flow  over  an  airfoil  and  can 
capture  the  rapid  evolution  of  shocks  in  this  compli¬ 
cated  flow.  A  maximum  of  224  interferograms  per  test 
condition  at  rates  of  up  to  40  KHz  can  be  recorded 
using  the  system.  To  date,  several  thousand  images 
have  been  obtained  at  rates  of  up  to  20  KHz.  A  com¬ 
prehensive  image  processing  software  package  has  also 
been  developed  for  processing  these  large  number  of 
interferograms.  The  software  has  reduced  the  time 
required  for  this  task  considerably.  It  includes  fea¬ 
tures  for  special  treatment  of  each  image  using  dig¬ 
ital  filtering,  contrast  enhancement,  local  magnifica¬ 
tion  for  high  fringe  density  regions,  accurate  fringe 
centerline  detection,  fringe  identification,  etc.  Inter¬ 
ferograms  containing  multiple  shocks  in  the  flow  have 
also  been  satisfactorily  analyzed.  The  effects  of  vari¬ 


ables  such  as  PDI  plate  optical  density,  pin-hole  size, 
etc.,  have  been  investigated  The  technique  is  robust, 
easy  to  set-up  and  offers  considerable  advantages  over 
other  interferometric  methods.  It  is  a  real-time,  self¬ 
aligning  technique  and  is  not  affected  by  extraneous 
influences  such  as  tunnel  vibrations. 
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Fig.  1.  Schematic  of  the  Compressible  Dynamic 
Stall  Facility. 
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Fig.  2.  Schematic  of  the  Point  Diffraction  Interferometry 
System. 
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Fig.  3,  Block  Diagram  of  Camera/Laser  Synchronization  for  the  High-Speed  PDl  System. 
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Fig.  4.  Flowchart  of  Events  Following  Handswitch  “FIRE”  for  the  High-Speed  Interferometry  System. 
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Fig.  5.  Timing  Sequence  for  the  High-Speed  Interferometry  System. 
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Fig.  8.  Representative  Interferograms  of  the  Row  Reid;  M  =  0.3,  k  =  0.05,  a  =  10.07  deg. 
(a)  Single  Exposure  Camera,  (b)  High-Speed  Camera. 


Fig.  9.  Comparison  of  the  Global  Pressure  Distributions  at  M  =  0.3,  k  =  0.05, 
a  =  13.99  deg  for  Different  Oscillation  Cycles. 


Rg.  10.  Effect  of  Pin-Hole  Size  on  PDI  Images;  M  =  0.4,  k  =  0,  a  =  8  deg. 
(a)  Small  Pin-Hole,  (b)  Large  Pin-Hole. 
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Fig.  12.  Automated  Image  Processing;  M  =  0.45,  k  =  0.05,  a  =  1 0.0  deg.  (a)  Digitally  Enlarged 
Region  of  PDI  Image,  (b)  Map  of  Dark  Fringes. 


